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A map from Id Quantum Field Theory to Quantum Chaos on a 2d Torus 
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Dynamics of a class of quantum field models on Id lattice in Heisenberg picture is mapped into a class 
of 'quantum chaotic' one-body systems on configurational 2d torus (or 2d lattice) in Schrodinger 
picture. Continuum field limit of the former corresponds to quasi-classical limit of the latter. 

PACS numbers: 05.45. +b, 05.30.Fk, 03.65.Fd, 72.10.Bg 



Understanding of long-time dynamics of interacting 
quantum many-body systems or quantized fields is a 
long standing open problem. In particular, one would 
like to understand the conditions for the emergence of 
quantum mixing (implying ergodicity), i.e. general de- 
cay of time autocorrelation functions. Inspired by rich 
quantum behaviour of non-integrable few body systems 
having (partially) chaotic classical limit ]jj few papers 
appeared recently concerning 'quantum chaos' in non- 
integrable many-body systems ^ |j. In Q dynamical 
phase transition from non-ergodic dynamics exhibiting 
ideal transport to mixing dynamics exhibiting normal 
transport has been demonstrated numerically in a non- 
integrable kicked t-V fermion model on Id lattice. 

Below we construct an exact linear mapping from 
dynamics of a certain large class of interacting infi- 
nite spin— i chains in Heisenberg picture to a class of 
non-linear One-Body Image dynamical systems (OBI) in 
Schrodinger picture which are realized either on config- 
urational 2d torus or on 2d lattice (tight-binding (TB) 
formulation). Further we will show by working out two 
examples: (i) how integrable behavior of the infinite XX 
spin chain in spatially modulated transversal magnetic 
field is connected to the integrability of OBI and to the 
Harper equation || , and (ii) demonstrate the phase tran- 
sition of the related non-integrable kicked XX chain from 
non-ergodic dynamics to ergodic and mixing dynamics 
corresponding to the stochastic transition from regular 
to chaotic motion of OBI in the classical limit. 

Let cr|, j e Z, s £ {x, y, z} denote a chain of indepen- 
dent Pauli spin— h variables. We start by generalizing 
the result of |J, namely we find that the operator space 
spanned by the following set of spatially modulated ob- 
servables {£/„(#), V n {d);n £ Z, £ [-tt, tt)} 



j=-oo 



j=-oo 



n > 1, 

—a] n = 0, 

o](o*)^o]_ n »<-l, 

a] (ojOn-i o- v j+n n>l, 
1 n = 0, 

-of (<T?)_ n _i cr 



j-n 



n < -1. 



where (oj)fc := ri;=i °f+i for fc > 1 and (cr|) := 1, is 
closed under the Lie bracket [A, B] = AB — BA and forms 
an infinitely dimensional Dynamical Lie Algebra (DLA) 



[U n (#),Ut(<p)] = 2iexp{%(W + ru P)at-n)V n -t(d + <p) 

- 2iexp(^(W + nip)s n -i) Vi- n (fl + <p), 
[U n (^),Vi(<P)] = 2£e X p(f(-W + n^)s i )C/„ +; (tf + ^) 

- 2iexp(%(W + mp)si) U n -i[& + <p), 

[V n {d),Vi{ip)] = (s n + sO{sin(i(M + n^)) x 

((s n -iSi + l)V n -l(0 + (p)~ {si-nSl + 1)V5_„(0 + <p)) 

+ 2sm(±(W-mp))V n+ i(d + tp)}. (1) 

where s n '■= — 1, 0, 1 for n <, —, > 0, resp., is a sign of in- 
teger n. Few notable members of DLA are: Ising or XX 
hamiltonian H 1 = JUi(0),H XX = J(Ui(0) + U-i(0)), 
spin interaction with modulated transversal magnetic 
field h z = hcos(ej) with period 27r/e lattice spacings 
H mh = \h(U (e) + U (-e)), spin current j s = Vi(0) + 
y_i(0), etc. Let us fix the fundamental field modulation 
e and introduce the following notation: 

(n, k) £ Z 2 , U± k = | (U n (ke) ± U_ n (-ke)) , 

V+ k = \(V n (ke) + V- n (-ke))s ni 

VTt = £ (V„(fce) - V_„(-fce)) , 



= U1 



iV, 



± 



,k> 



(y,x)£T 2 , W ± {y,x) 



2vr 



E i(ny+kx)yrr± 



We may also consider e as a lattice spacing and treat 
W^ k as a set of spatially 27r-periodic fields. We will 
assume that the modulation is incommensurable with 
the lattice spacing, i.e. that e/27r is irrational, other- 
wise obsrvables W„ fc are periodic w.r.t. index k. DLA 
becomes a Hilbert space when we introduce an infi- 
nite temperature (grand) canonical scalar product Q 
(A\B) := limi^oo j-2~ L tr A^B, where L is a diverging 
length of the spin-chain. Let the two linear subspaces 
spanned by W° k (or W a (y,x)) for a £ {+, — } be de- 
noted by SCTlo-. The spaces and SDT_ are orthogonal 
and observables W^ k (n,k £ Z) and W~ k (for n > 1) 
form orthonormal bases in each of them, since one can 
show 

(Kj W Xl) = W*,«> {W^AW-J = 0, 

(W-JW- ,) = (5 n , m - S n ,. m )S ktl . (2) 



The full set {W n k } is over-complete, since W_ n 



-W. 



n,k ' 



while the subspace 971+ = 971+ is self-adjoint, 

since \ — W^ n _ k . One can write analogous rela- 
tions in terms of continuous variables (y,x). We have 
DLA = Ml © 971+ © 971- . Note that the adjoint map 
(ad A)B = [A, B] generates the Heisenberg motion on 
DLA, exp(it&dA)B = e itA Be~ ltA . In particular, the mo- 
tion generated by k has a beautiful structure. Let us 
write the self-adjoint Hamiltonian in a general form as 



n,k 



(3) 



using two sets of possibly time dependent real coef- 
ficients g n ,k = g n ,k(t),7n,k = ln,k{t). Tedious but 
straightforward calculation, using algebra (|j), gives the 
action of ad_ff on two continuous sets of observables 
W + (y,x),W~(y,x), (y,x) e T 2 which can be written 
in terms of two non-local 'Schrodinger operators' fl^ 



TV 



(adif) W ± (y,x)=-tH ± W ± (y,x), (4) 

^ %n,fe (sin (np x - kp y ) sin (foe + ny - 7„ )fe ) 
n,k - sm(np x + kp y )sm(kx-riy-7 nt k)), 

hg 7hk (cos(np x - kp y ) cos(kx+ny-j n ,k) 
n > k + cos(np x + kp y ) cos(kx-ny-j n ,k)) , (5) 



where p^j, = —i%d/d x>y are momentum operators conju- 
gate to £,y with an 'effective Planck constant' 



Since Heisenberg dynamics generated by H is closed on 
SDTcr, (adi/)97t (T C 97l CT , one may write a general time- 
evolving operator A(t) e s XSl a in terms of a complex- 
valued 'Schrodinger wave function', in either 'momen- 
tum' ^ A k (t) or 'position' ty A (y,x;t) representation 



A{t) 



n . k 



KM W n,k = / dydxy A (y,x;tyW°(y,x). 



T- 



By means of eq. (|j) and the fact that H a is Hermitian 
on L 2 (T 2 ) (which can be checked directly using the ex- 
pressions (p)) one can easily show that the Heisenberg 
evolution of the observable A(t), (d/dt)A(t) = i(&dH)A, 
is fully equivalent to the Scrodinger equation 



ih-* A (y,x;t) =H°* A (y,x;t). 



(7) 



governing time evolution of one particle on a torus T 2 
(OBI). The bilinear map (H, A(t)) «-> (H a ,^ A (y,x;t)) 
is a central result of this Letter. To conclude a general 
exposition we make few remarks: (i) A non-trivial 'classi- 
cal limit' h — > of OBI exists, being equivalent (^) to the 
continuum field limit of the quantum spin chain model 



e — ► 0, if hg n j (and not g n ^i alone) are kept constant and 
finite, (ii) The operators Tt + and H~ (|J) commute 

[H+,H-} = 0, (8) 

and the Poisson bracket of the corresponding classical 
counterparts vanishes, (iii) As a consequence of the pre- 
vious remark we find that OBI ^ (and its classical limit) 
is integrable, 7i _CT being the second integral of motion, 
provided the original spin- field Hamiltonian H (g) or OBI 
Hamiltonian H a is autonomous, i.e. (d/dt)H = 0. How- 
ever, one has a possibility of chaotic motion in classical 
limit and emergence of 'quantum chaos' when the prob- 
lem is explicitly time-dependent, say that coefficients are 
periodic functions, g n ,k(t + 1) = 9n,k(t)- In such case one 
integrates the evolution over one period of time and de- 
fines the unitary Floquet maps U = Texp(— i J dtH(t)), 
U a = 7~cxp(— i L dfH" (t)) . (iv) Temporal correlation 
functions of the quantum field problem are mapped (us- 
ing eqs. (g)) onto transition amplitudes of OBI 



(A(t)\B(t')) 



(y A (t)\V B (t')) 

{^ A {t)\Vy^ B (t')) 



(9) 



where T' v ^{y 1 x) = ^(y,x) — y, x). Therefore, er- 
godic properties of many-body dynamics on DLA are 
determined by the spectral properties of OBI: (a) Spin 
chain is quantum mixing in 97to-, lim t _ >00 (A(t)|_B) = 
0, A, B £ Wl a , iff the spectrum of OBI Hamiltonian 
H' 7 (or of OBI Floquet map U") does not have (non- 
trivial) point component, (b) Spin chain is quantum er- 



( 6 ) godic in 97! CT , lim T ^oo T^ 1 J Q T dt(A(t)\B) = 0, iff (or 
1) is not in the non-trivial point spectrum of TV 7 (or 
If 7 ). In autonomous case, dH/dt = 0, the Hamiltonian 
H and the trivial zero-frequency eigenstate, 7i +1 i' H = 



. m lH,H] = 



g-n,ke 



' +9- 



k 8 



h(9n,ke~ 



+ g n - k e l 



'•) should be excluded from 



971+ and L 2 (T 2 ), respectively, i.e. (A\H) = (B\H) = 0. 

We apply the above results to work out two inter- 
esting examples. Example I: XX spin chain in spa- 
tially modulated quasi-periodic transversal magnetic field 
h 3 = (0,0,/icos(ej)) (XXmh) 

H = H xx + H mh = JU+ + ^h(U+ A + U+_ x ). (10) 

Here the Heisenberg dynamics on DLA is governed by 
the following commuting one-body problems 

"H + = a $mp x siny — psin.p y since, (11) 
HT = a cosp x cos y + (3 cosp y cos x, 

where a = 2eJ = 4HJ, (3 = 2eh = 4Hh. This models are 
directly related to the electron motion on 2d rectangular 
a x 6 lattice in a uniform perpendicular magnetic field 
h' within the TB approximation |g,^|. In the symmetric 
gauge A — ~h'(—y, x, 0) the TB problem with the band 
energy £{K) = § cos(a_ftTi) + % cos{bK2) reads 



+ I ( e -H«"* nifc+1 + e^Xk-i) (12) 

where e = e a abh' / c Q ?iphys |j is here the dimensionless 
magnetic flux thru lattice cell. We note that discrete 
indices (n, k) £ Z 2 now label the position lattice (na, kb) 
while continuous indices (y, x) G T 2 are the conjugate 
quasi-momenta. OBI Hamiltonians TL^ can be written in 
terms of TL and its time-reversalTL* = H.\h'^-h> , namely 
W ± = HtT-L*, and hence [H,H*] = 0. Using a different, 
Landau gauge A = h'(0,x,0) the TB problem (|l2) can 
be re-written in terms of Id Harper equation [|5| 



ka(u n +i + u n _i) + /3cos(ne — = J5tt r , 



(13) 



Let us assume for the moment that a < [3. Then 
Uni'&'iE) — u n is a unique exponentially localized eigen- 
function (EF) of eq. ( filf ) which has a dense pure point 
spectrum, and ^ n ^{d]E) = exp — |en)fc) -E) 
is a degenerate dense set of EFs of TB problem (12), 
E) — Ety n .E), for a dense set of parame- 
ters i? H . Though 7i and 7i* should have a common set 
of EFs, * n ,fc(^; E) is not an EF of ft*, neither it is in L 2 
since it is extended in variable k. We search for such EF 
with an ansatz 3> nj fc($; E, E') = J^j v j^n,k{$+£]', E) and 
require ft*$ n , fe (i?; E, E') = E'$ n)fe (i?; E, E') yielding the 
Harper equation ( fi"3| ) for coefficients w„ = u n (d;E'). 
Thus we obtain a common set of EFs of ft and ft* in 
terms of a 'convolution' of two Id Harper functions 



i(i?-d-i£ra)fc 



*„ lk (0; E, E')=J2u l+n {&, E)u0- E')e 



which is also a common set of EFs of Ti.^ ( |Tl| ) 
H^nA^ E , E ') = ( E T S')*»,fc(^ 



(14) 



(15) 



The property + e;E,E') = <S> n , k ($;E,E') sug- 

gests independence of EF on parameter ■& provided e/27r 
is irrational. If a > (3 localized EFs can be constructed 
analoglously by 'duality transformation' n <-> k,y <-> x. 
Thus, we found that OBI 7i ± have a dense pure point 
spectrum for a ^ (3, hence XXmh is non-mixing, non- 
ergodic and even completely integrable, namely 'zero- 
energy' eigenstates of Ji a are the images of a complete set 
of conserved charges, Q a {E) = J2 n ,k®n,k( E i aE ) W Z,k> 
[H,Q a (E)] = 0. However, time autocorrelation func- 
tion (j4(t)|A) of a certain observable A may still de- 
cay to zero provided the image function ^ A is or- 
thogonal to all (localized) EFs ([l4j) of H ± . Interest- 
ingly, this happens with the spin-current j s — W^~ — 
W± lfi if a < (3 (i.e. J < h), since EFs (0) have 
the following properties <^ n ^{E,E')* = $„.fc(E ', E) = 
and $ n , (E,E')* = $ n , (E,E') (we put 
:= 0) implying $„, (E, E') = ®- nfi (E, E'). So we have 
{W-\$(E,E')) = 0. This proves (j„(t -> oo)|j s ) -» 



and non-ballistic spin-transport (vanishing spin stiffness 
D s := lim T ^ 00 (l/T) £ dt(j s (t)\j s )) for J < h, while for 
J > h we find in general ballistic transport (D s > 0) 
since no similar symmetry exists for the other index fc. 

Example II: kicked XXmh model (kXXmh) with time- 
dependent Hamiltonian 



H(t) = JU+ + ^h(U+ 1 + U+_,) £ 6(t 



(16) 



One-period propagator from just after the kick 

U = exp (-qh(U+ tl + U+_ 1 )) exp (-iJU+ Q ) (17) 

is equivalent to Floquet quantum maps of two kicked OBI 

U + = exp(^ sinpy sinx) exp sinp^. siny) , (18) 
U~ = cxp(^p cospy cosx) exp cosp x cosy) . 

In the following we will consider only the map U + since 
the space 971+ contain physically more interesting observ- 
ables. The Floquet evolution * A (m) = U+ m ^ A {0) yield- 
ing Heisenberg evolution of observables A(m) G is 
in the 'classical' limit equivalent to a volume-preserving 
(2 x 2)d map on T 2 x M 2 

x' = x + a cos Pa sin y, p' = p y — a sinp^ cos y, (19) 
y' = y — (3 cos p' y sin x' , p' x = p x + /3 sinp' y cos x , 

which is non-integrable and (almost) fully chaotic for suf- 
ficiently large kick parameters, a, (3 3> 1. Interesting 
question is now if and when the dynamics of kXXmh is 
quantum mixing and how it corresponds to dynamics of 
the 'classical map' ( |T9| ) as h = \e — > 0. This problem has 
been approached numerically by iterating the one-body 
Floquet map IA + on a finite (truncated) momentum space 
(n, k) e {-N/2 . ..N/2} 2 . The position states are then 
discretized as Xj — sj,yj — sj.s — 2tt/N. The trun- 
cated Floquet map U + can be efficiently implemented by 
means of Fast Fourier Transformation (FFT), namely if 
F is Id FFT on N sites then N 2 x N 2 Floquet matrix is 
de-composed as 

{F- 1 ® l)(diag(7 n , fe )(F ® F- 1 )(diagZ?„, fe )(l ® F), 
with diagonal matrices C n> u = exp (i((3/h) sin sn sin hk) 
and D n ,k — exp(—i(a/h)smhnsmsk), requiring ~ 
AN 2 log 2 N computer operations per time step. In order 
to avoid recurrences of quantum probability due to finite- 
ness of momentum space we use an absorbing boundary 
in momentum space, namely after each iteration of the 
truncated Floquet map wc multiply the wave- function by 
a box-window, *„,fc(m) -> 9(N/2 - a/h - \n\)0(N/2 - 
(3/h — |fc|)^' n; fe(m), Convergence to true dynamics on a 
torus has been checked by comparing results for different 
truncations, say ./V and N/2 (we went up to N = 2 14 ). 

In Fig.l we show numerical results for the auto- 
correlation function of the spin current C(m) = 
(J s (m) \j s (Q)) while similar, compatible results have been 



obtained for the time-correlations of other observables. 
(i) For sufficiently large kick parameters a, (3 the classi- 
cal map (19) is strongly chaotic and mixing exhibiting 
normal diffusion in momentum plane (j> x ,p y ). However, 
kXXmh is not exactly mixing for any finite h: \C(m)\ is 
rapidly (possibly exponentially) decreasing down to some 
value C* — \C(m)\ where it saturates. When we de- 
crease h, C* decreases proportionally C* tx h, and so 
in the 'quasi-classical'/continuum limit H = — > the 
point spectrum of U + vanishes and kXXmh approaches 
mixing behaviour in accordance with the map (|l9[). (ii) 
For smaller but still finite values of a,j3 the classical 
map enters into the regime of KAM quasi-integrability 
with invariant tori suppressing the diffusion of momenta 
(p x ,p v ). Correspondingly kXXmh is non- mixing and 
C* ~ 1 for any value of h. In this regime, C{m) is very 
weakly h— dependent. In both regimes, (i) and (ii), the 
square widths of the 'wavepackets' (^ 3a (m)\p%. Jv^ 3 ' (to)) 
have been found to be uniformly increasing in time and 
limited only by the size of the truncated momentum space 
N. This rules out the possibility of quantum localization 
and existence of pure point spectrum, and indicates co- 
existence of point and continuous spectrum for any finite 
h (and finite a, a situation similar to (possibly re- 
lated) Id kicked Harper model ||. In the limit a, J3 — > 0, 
the continuous spectral component vanishes and we re- 
cover integrable XXmh model with pure point spectrum 
as discussed above. The quantum correlation function 
C(m) seem to follow the quasi-classical propagator only 
up to logarithmically short time, namely we found em- 
pirically that deviation (when it is small) increases expo- 
nentially \C(m)-C n ^o(m)\ » 0.022ft 2 e Am with A « 0.59 
for a = 3,(3 = 0.75 and A w 1.1 for a = 6,0= 1.5. 

Conclusions. In a specific ood class of (Pauli spin, or 
spinless fermion) quantum field models in Id, the Heisen- 
berg time evolution in two disjoint ood linear subspaces 
of essential field observables has been shown to be for- 
mally equivalent to the Schrodinger dynamics of a class 
of one-body image problems on a 2d torus (or 2d lat- 
tice) . Autonomous models of this class were found to be 
completely integrable, pointing out a novel class of inte- 
grable one-body problems (|^,^). For example, dynamics 
of XX chain in a static quasiperiodic transversal field has 
been solved in terms of Harper equation JTcj ] . However, 
time-dependent (e.g. periodically kicked) models of our 
class behave in a non-integrable fashion being mapped 
onto one-body problems with chaotic classical limit. It 
seems that spatial modulation is crucial to break integra- 
bility since spin chain kicked with homogeneous transver- 
sal field remains completely integrable as found in || . In 
the contunuum field limit our kicked spin chain model 
(kXXmh) has been demonstrated to undergo a (phase) 
transition from mixing to non-mixing dynamics (similar 
to a transition found in |Q), as its one-body counter- 
part in the classical limit undergoes a stochastic transi- 



tion from chaotic to quasi-regular motion. This is an in- 
teresting link between quantum field theory and chaotic 
dynamics and should inspire future research in this di- 
rection. Such approach to long-time dynamics of certain 
(non-integrable) quantum many-body systems, since it 
makes time evolution formally equivalent to 'quantum 
chaos' in few degrees of freedom, overcomes the tradi- 
tional problems due to huge Fock space in thermody- 
namic limit. Financial support by the Ministry of Science 
and Technology of R. Slovenia is acknowledged. 
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FIG. 1. In |C(m) in kXXmh for several indicated values of 
parameters a and e while j3 — a/4. Note the transition to 
mixing dynamics as e is decreased (heavy-light full curves) 
for 'chaotic case' a — 6, and stable non-mixing behaviour 
for 'quasi-regular' case a = 3 (dotted curve, curves for other 
(small) values of e are almost indistinguishable). Broken lines 
at In | e/3 indicate the scaling C* oc !i = he. In two insets 
we show two orbits in momentum plane (chaotiCjdiffusive for 
a = 6 and quasi-regular for a = 3) of the map (Um of length 
3000 starting at xq = 0.2, yo = 0,p x o = 0,p v o = 0.005. 



